Uniaxial-strain detwinning of CaFe 2 As 2 and BaFe 2 As 2 : optical and transport study 
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The parent compounds of iron-arsenide superconductors, y4Fe2As2 (A=C&, Sr, Ba), undergo a 
tetragonal to orthorhombic structural transition at a temperature Tto in the range 135 to 205 K 
depending on the alkaline earth element. Below Tto the free standing crystals split into equally 
populated structural domains, which mask intrinsic, in-plane, anisotropic properties of the materials. 
Here we demonstrate a way of mechanically detwinning CaFe2As2 and BaFe2As2. The detwinning is 
nearly complete, as demonstrated by polarized light imaging and synchrotron X-ray measurements, 
and reversible, with twin pattern restored after strain release. Electrical resistivity measurements in 
the twinned and detwinned states show that resistivity, p, decreases along the orthorhombic a -axis 
but increases along the orthorhombic & -axis in both compounds. Immediately below Tto the ratio 
Pbo/Pao = 1-2 and 1.5 for Ca and Ba compounds, respectively. Contrary to CaFe2As2, BaFe2As2 
reveals an anisotropy in the nominally tetragonal phase, suggesting that either fluctuations play a 
larger role above Tto in BaFe2As2 than in CaFe2As2, or that there is a higher temperature crossover 
or phase transition. 

PACS numbers: 74.70. Dd,72.15.-v,68.37.-d,61.05.cp 



INTRODUCTION 



The parent compounds of high-transition temperature, 
T c , iron arsenide superconductors, i?FeAsO (R stands 
for rare earth) [l[ and ^4Fe2As2 (A = Ca, Sr, Ba, we 
denote compounds as A122 in the following) 0, 0] un- 
dergo a tetragonal to orthorhombic structural transition 
upon cooling below a transition temperature, Ttoj a c- 
companied or followed by an antiferromagnetic ordering 
at Tm 0. In free standing single crystals, this transition 
leads to the formation of twin domains of four types 
In Ba(Fci_ 2; Co 2 ;)2As2, Tto decreases with increasing x 
and reaches T c for x «0.063 For a range of dop- 

ings twin domains coexist with superconductivity and 
strongly affect it [1, ■ While the structural anisotropy 
(the ratio of the in-plane lattice parameters, a/b) is very 
small, the electronic structure in the AFM phase is highly 
anisot ropic, as manifested in the first principles calcula- 
tions [lll - ll4l | and in the experiment [la]. In particular, 
the Fermi surface lacks even approximate 4-fold symme- 
try. At the same time, despite a dramatic anisotropy in 
the Fermi surface geometry, calculations predict a rather 
small transport anisotropy (see discussion below). In- 
trinsic strong anisotropy is a key component in a popular 
theoretical picture of "nematic" antiferromagnetic state 
(see review Refflrjl for a discussion) , which is believed to 
be incipient even outside the long-range ordered antifer- 
romagnetic phase, and probably plays an important role 
for superconductivity [F6|-ll9(. Yet, experimentally it has 
not been accessible till now, because of the twinned do- 
main structure masking the internal anisotropy. There- 
fore, it is important to reliably obtain and characterize 



textured, or even better, single domain samples of theses 
materials. 

The fact that Tto is below room temperature, makes 
the situation very different from the classical case of the 
cupratc YBa 2 Cu307_,5, where single orientations of the 
orthorhombic domains can be stable during sample han- 
dling and over the entire experimental temperature range 
studied [2(j. Therefore different approaches are needed 
to reach a single domain state. The first attempt to study 
intrinsic in-plane anisotropy in the iron arsenides was 
utilizing a strong magnetic field to align domains [2l| . 
similar to the effect used for turning La2Cu04 into a 
single-domain state [22|. The application of a magnetic 



field of 14 T led to a change of the domain patterns in 
Ba(Fe 1 _ a; Co a ;)2As2 (£=0.016 and 0.025) and in sample 
with £=0.025 lead to an increase of one of the popula- 
tions by 7 %, (from 54 to 61%.) This resulted in a non- 
sinusoidal angular dependence of the magnetoresistance 
with a magnitude of several percent, and, because of a 
small population imbalance, was taken as a signature of 
large in-plane resistivity anisotropy. 

Here we report nearly complete reversible mechani- 
cal detwinning in the parent compounds CaFc2As2 and 
BaFe2As2- We used uniaxial strain to drive a region of 
the crystal into a single domain state and study the tem- 
perature dependence of its electrical resistivity. In both 
CaFe2As2 and BaFe2As2 we find a clear effect of detwin- 
ning on the resistivity below Tto> with a resistivity de- 
crease along the orthorhombic a Q -direction and an in- 
crease along the b a direction. Here and below we use the 
notations where a a > b 0l so that the ferromagnetic chains 
run along b . That is to say, the observed anisotropy is 
not only much smaller than one may anticipate from 
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a naive orbital-ordering picture, but also the sign of the 
anisotropy is opposite to what one could expect: the elec- 
trons move easier along the antiferromagnetic direction 
than along the ferromagnetic one. Interestingly, both 
these results also appear in the first principle calcula- 
tions. 

The ratio Pbo/Pao is maximal right below TtOj an< i 
is larger in Bal22 than in Cal22 (1.5 vs 1.2). Despite 
the fact that the magnetic moment and the degree of or- 
thorhombicity are increasing upon cooling, the transport 
anisotropy is decreasing. We discuss possible explanation 
of this anomalous behavior below. 

Finally, the resistivity of Bal22, but not of Cal22 
shows measurable anisotropy for at least 30 K above Xto, 
reflecting nematic fluctuations above Xro [IH in the for- 
mer compound, consistent with the fact that the phase 
transition in Cal22 is strongly 1st order, and in Bal22 
2nd order or close to that. 



METHODS 
Experimental 

White-light, optical images were taken at tempera- 
tures down to 5 K using a polarization microscope Leica 
DMLM equipped with the flow- type 4 He cryostat, as de- 
scribed in detail in Ref. High-resolution static images 
were recorded. The spatial resolution of the technique is 
about 1 fxm. Single crystals of Cal22 were grown from 
Sn flux as described elsewhere. 3, 25| Single crystals of 
Bal22 were grown from FeAs flux [7[. 

As schematically shown in Fig. [IJ panels (a) and (b), 
four types of domains are formed in the orthorhombic 
phase due to orientational degeneracy of the direction 
of orthorhombic distortion Since domains 01 and 
02 (and similarly 03 and 04) share a common plane 
corresponding to the tetragonal (100) [or, equivalently 
(010) plane], their formation does not require lattice de- 
formation and they easily form pairs. This leads to a 
slight deviation of the orthorhombic axes in the pair from 
perfectly mutually orthogonal orientation, by an angle 
±a = 7r/2 — 2 arctan (b/a). For resistivity anisotropy 
study, it is important that domains with two orthogonal 
orientations of orthorhombic a Q axis are pairwise inter- 
mixed, averaging the in-plane anisotropy in the twinned 
state. 

Straining the crystal along the tetragonal [110] direc- 
tion (that become orthorhombic a or b Q axes in the or- 
thorhombic phase) removes the degeneracy and leads to 
preferential orientation of domains with the a D axis along 
the strain. Technically, this requires application of strain 
at 45° to the natural tetragonal [100] and [010] facets of 
the crystal [2f|. Due to both very clean natural growth 
faces and difficulty to form a good cleave surface, in our 
initial attempt to dctwin samples of CaFe2As2 we used 



as-grown samples. Mechanical strain was applied to the 
sample through thick silver wires (125 /xm thick), sol- 
dered with tin-based alloy (2?| to the corners of the sam- 
ple. Use of silver wires to transmit the strain allows for a 
gentler sample deformation, which is very important for 
a material as soft as CaFe2As2- It also allows electrical 
contact to the sample for electrical resistivity measure- 
ments. The strain was applied by mechanically deform- 
ing a brass horseshoe with a stainless screw, Fig.QJc), at 
room temperature. The wires were soldered to the pre- 
strained horseshoe, so that both strain and a little stress 
(due to wire stiffness) could be applied. 

The domain population of the BaFe2As2 sample has 
been analyzed by high-energy x-ray diffraction. Entire 
reciprocal planes were recorded using the method de- 
scribed in detail in Refl28l which has been successfully 
applied recently to study the domain structure in pnic- 
tides 0, IH. The absorption length of the high energy 
(99.3 keV) x-rays from the synchrotron source (beamlinc 
6ID-D in the MUCAT sector at the Advanced Photon 
Source, Argonnc) was about 1.5 mm. This allowed full 
penetration through the roughly 0.3 mm thick sample 
mounted with its c direction parallel to the incident x- 
ray beam, which was reduced to 0.2x0.2 mm 2 size by 
a slit system. Therefore, as result each single measure- 
ment averages over the entire sample volume selected by 
the beam dimension in the (ab) plane and its projec- 
tion through the sample along the c direction. The di- 
rect beam was blocked by a beam stop behind the sam- 
ple. Two-dimensional scattering patterns were measured 
by a MAR345 image-plate positioned 1730 mm behind 
the sample. During the recording, the sample was tilted 
through two independent angles, [i and 77, perpendicular 
to the incident x-ray beam by 3.2 deg. 



Theoretical 

First principles calculations were performed using the 
standard LAPW method as implemented in the WIEN2k 
package (29| . The experimental lattice and internal pa- 
rameters were used for both compounds. In order to 
estimate the transport anisotropy within the constant re- 
laxation time Boltzmann approximation, we have calcu- 
lated the anisotropic plasma frequency, ujp a = 4ir (Nv^) , 
where TV is the density of states per unit volume, a = x, y, 
(with x along a a direction), v a is a projection of the Fermi 
velocity and the brackets denote averaging over the Fermi 
surface. Up to 70000 k-points in the Brillouin zone have 
been used to perform this averaging. 

To test the sensitivity of the results, we have performed 
calculations in both Local Density approximation (LDA), 
with smaller magnetic moments and larger Fermi sur- 
faces, and Generalized Gradient approximation (GGA), 
with larger magnetic moments and smaller Fermi sur- 
faces. 
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RESULTS 
Detwinning 

The series of pictures in Fig. [1] shows polarized opti- 
cal images of the Cal22 sample, taken at 5 K, in the 
area between the contacts without external strain (d), 
at progressively increasing strain, (c) through (g), and 
after strain release (h). On the application of strain 
the domain pattern in the central area between contacts 
changes rapidly and eventually a domain free region is 
formed. However, a pattern of randomly oriented do- 
mains is still seen away from the central part, reflecting 
the distribution of strain over the sample area and its con- 
centration close to the soldered contacts. Strain release 
restores the domain pattern, i.e. detwinning is reversible. 

Encouraged by this initial success, we introduced two 
modifications into the sample detwinning process. First, 
samples were cut into strips to allow a homogeneous dis- 
tribution of strain between contacts. The image of the 
sample after cutting, Fig. HJa) , reveals a pattern of ran- 
dom domains, as is typical for an unstrained sample. Sec- 
ond, we attached contacts in a conventional four-probe 
configuration. The thick wire contacts were made at the 
ends of the sample and served as current contacts. The 
potential leads were made from a much thinner, and thus 
less rigid, 50 /im silver wire long enough to avoid gener- 
ation of additional strain. The potential contacts were 
mounted on the rear surface of the crystal, so that we 
could monitor the detwinning in a polarized optical mi- 
croscope. Fig. &b) shows an optical image of a sample 
strained through current contacts. Although there is a 
clear change in the domain distribution upon the applica- 
tion of strain, showing that detwinning affects the whole 
sample thickness (30l |. two big unstrained areas can be 
seen in the central region, which corresponds to the ar- 
eas above the potential contacts on the rear surface. This 
observation clearly shows that the contacts affect strain 
distribution in the sample and do not allow homogeneous 
sample deformation. 

In an attempt to further improve detwinning for re- 
sistivity measurements, we started sample deformation 
through potential contacts. In this case only a central 
part of the sample, between the potential contacts, is de- 
formed. Simultaneously, the parts of the sample between 
the thin- wire current contacts and the thick-wire poten- 
tial contacts remain strain-free. Imaging of the sample at 
5 K during different stages of the sample preparation and 
straining is shown in Fig. [3] As can be seen in a series of 
images, the central part of the sample is predominantly 
free of domains, while a clear pattern of domains can be 
seen in the unstrained part of the sample. The pattern of 
domains disappears abruptly at Tto, the strained part 
abruptly changes color at Ttoj signaling the disappear- 
ance of the orthorhombic distortion. Formation of the 




FIG. 1: (Color online) (a,b) Schematics of domain formation 
during tetra-ortho transition. Orthorhombic phase is real- 
ized through displacement of atoms along the diagonal of the 
tetragonal lattice with doubling the unit cell volume and 45° 
rotation of the crystallographic axes [25|. The domains are 
distinguished by the direction of the distortion, a , shown for 
four domains with the arrows, (c) A single crystal sample 
of CaFe2Fe2, mounted on a horseshoe, (d-h) Polarized light 
image of the crystal without external strain (d), strained to a 
progressively higher values (e-g), and after strain release (h). 



orthorhombic phase is homogeneous through most of the 
sample, except for two small patches close to the con- 
tacts. 
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FIG. 2: (Color online) A polarized light image (5 K, T « 
Tro) of a single-crystal strip of CaFe2 As2 cut along the tetrag- 
onal [110] direction before contacts were made (top). The 
sample was mounted in a usual four probe contact configu- 
ration with potential contacts located on the sample surface 
opposite to the surface in the image. Strain was applied to the 
sample at room temperature through current contacts and the 
image was taken at 5 K on a clear surface under the poten- 
tial contacts. A big single-domain area can be seen between 
potential contacts, however, large heavily twinned areas are 
found above the contacts (circles in the image), revealing no- 
table local reduction of strain by surface tension at the contact 
points. 



Resistivity of CaFe2As2 

The left panel of Fig. |4] shows the temperature- 
dependent electrical resistance, measured on sample A 
during two successive thermal runs (cooling from 300 K 
to 5 K and warming up to 300 K) in the unstressed (blue 
curve) and stressed (red curve) states. It is clear that 
the strain itself is not large enough to cause a change in 
the value of the sample resistance at room temperature, 
beyond a small (~ 1%) systematic resistance increase on 
strain change (both increase and decrease) due to a fa- 
tigue. The data above Tro stay unchanged, below a clear 
change in resistance is observed. In particular, the up- 
ward jump in the resistivity below Tro is notably dimin- 
ished in the detwinned state. The data in the unstrained 
state compare well with the temperature-dependent resis- 
tivity in the standard resistance measurements along the 
[100] direction @, 0, ED, El . The inset in the left panel 
shows the same data, normalized by the room temper- 
ature value, in comparison with the data for sample B, 
measured in the same contact configuration. The data 




FIG. 3: (Color online) (a) A polarized light image (5 K, 
T <C Tro) of a single crystalline strip sample of CaFe2As2 
with soldered contacts before being mounted on a horseshoe. 
On application of strain in a horizontal direction through the 
potential contacts, a strain-free area at the left end of the 
sample remains heavily twinned (b), while the central area 
turns into a nearly single domain state, with small twinned 
areas at the bottom of the image (c). After three thermal 
cycles from the 300 K to 5 K range, strain in the sample is 
partially released (d), however, large single-domain areas re- 
main, demonstrating the reproducibility of the strain-induced 
detwinning. 



agree very well in the magnitude of the effect. It can 
be clearly seen that the position of the transition is not 
shifted with strain, the transition remains sharp, and the 
width of the hysteresis between the cooling and warm- 
ing cycles remains the same. This clearly shows that the 
transition itself is insensitive to the applied strain, and, 
in view of the strong sensitivity of the transition tern- 
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FIG. 4: (Color online) Left panel. Temperature-dependent 
electrical resistance measured on a sample A of CaFe2As2 
in twinned (Rt, blue curve) and detwinned (R ao , red curve) 
states. The third curve is obtained as hypothetical Rt o — 
2R t — R a<1 , impossible to measure in the same contact ar- 
rangement, see text for details. Inset shows zoom of the 
transition region for samples A and B, revealing diminish- 
ing amplitude of the resistance jump in p ao (T) with identical 
transition sharpness and position. Right panel shows resis- 
tance of sample C, in which resistivity was measured along 
b-axis in transverse to strain contact configuration (see Fig. O 
in stress-free and stressed (partially detwinned) states. Inset 
shows zoom of the transition area. This data directly show 
an increase of the resistance jump in pb a . 



To get an independent assessment of the Pbo{T) we 
have mounted the sample C in such a way that the strain 
be generated in the direction perpendicular to the cur- 
rent. This was achieved by soldering the current con- 
tacts along the entire side of the rectangular sample and 
straining the sample and the current contact wires on 
a horseshoe. The ends of the sample were covered with 
epoxy glue and glued to two pieces of fiberglass to provide 
homogeneous deformation (sec Fig. [5]). The area of the 
potential contacts was minimized to diminish their effect 
on the measurements. Sample imaging showed clear de- 
twinning, see Fig.[5j however not in the entire sample. In 
the right panel of Fig. 2] we show the resistivity obtained 
in the twinned and partially detwinned states of sample 
C. The resistivity indeed shows an increase right below 
Tto At the same time, the temperature-dependent resis- 
tance in the tetragonal phase is unaffected by the strain 
all the way down to Tto, indicating that there is no resid- 
ual nematic ordering above Tto, which is, in fact consis- 
tent with the fact that the phase transition in Cal22 is 
strongly 1st order [l||. Note that since all measurements 
are taken on the same sample with the same contacts, 
they are free of any possible geometric uncertainties. 



Resistivity of BaFe2As2 



pcrature to uniaxial stress 3l| , the applied strain is very 
small and homogeneous. 

To check if the effects of thermal cycling are impor- 
tant, we took a polarized light image of the sample A 
after a resistivity run in the detwinned state. The op- 
tical image of the sample in a third successive thermal 
run, following the optical run of Fig. EJc) and resistivity 
measurements in Fig. 01 is shown in a panel (d) in Fig. EI 
As can be seen, the strain is slightly released after ther- 
mal cycling, however most of the sample remains domain 
free. After this second imaging, the strain in the sam- 
ple was released, and temperature-dependent resistivity 
measurements were undertaken in a strain-released state, 
main panel in the left panel of Fig. 2) These revealed 
a resistivity curve very close to the strain-free samples, 
restoring the magnitude of the resistivity up-jump to its 
initial value. 

Since the contact configuration in Fig. [3] does not al- 
low for a measurement of pbo, we inferred its evolution 
by assuming that the resistivity in the twinned state is 
an average of p ao and pb - This is strictly true if the 
population of the domains is random. The calculated 
Rbo(T) = 2R t — R ao is shown in the left panel. The com- 
parison of Pbo{T) and p a o(T) gives a ratio of 1.2 at the 
transition. With cooling, the ratio decreases to 1.05-1.1 
(due to the small absolute value, the effect of fatigue can- 
not be neglected, and it is hard to place a more precise 
value on the low-temperature anisotropy). 



Single crystals of Bal22 were detwinned in the same 
fashion as Cal22 above, with optical imaging confirming 
a homogeneous detwinned state for resistivity measure- 
ments along the a D -axis and a predominantly detwinned 
state for resistivity measurements along the 6 Q -axis. Be- 
cause of small image contrast, samples were additionally 
studied with high-energy synchrotron x-ray as described 
in the next section. 

In Fig. H3 we show the temperature-dependent electri- 
cal resistivity of Bal22. The left panel shows the re- 
sistance measured along the strain, R ao (sample D), the 
inset shows the transition area with the data from sample 
D as well as from another sample, E. The partial super- 
conductivity seen in these samples is presumably due to 
surface strain [33| associated with cleaving and shaping 
the sample and is not focus of this study. As for the 
case of CaFe2As2, we can infer the resistivity along the 
b a axis from Rb = 2R t — R ao ', the inferred Rbo is also 
shown. The right panel shows resistivity measurements 
with current transverse to the strain in the twinned and 
partially detwinned states on yet another sample F. 

Below Tto the data from Bal22 is qualitatively simi- 
lar to that of Cal22; both show anisotropies right below 
T TO , of 1.2 (Cal22) and 1.5 (Bal22). In both materials 
the anisotropy decreases with cooling below Tto, that is 
to say, is anticorrelated with the degree of orthorhom- 
bicity and the long-range magnetic moment. The abso- 
lute change of the anisotropy is similar, 10% (from 50% 
to 40% in Bal22, from 20% to 10% in Cal22). Note 
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FIG. 6: (Color online) Left panel. Temperature-dependent 
electrical resistance measured on a sample D of BaFe2As2 
in twinned (Rt, black curve) and detwinned R ao , red curve) 
states. The third curve is obtained as hypothetical Rbo = 
2R t — Ra , impossible to measure in the same contact arrange- 
ment. Inset shows zoom of the transition region for samples 
D and E. Right panel shows resistance for sample F, in which 
resistivity was measured along fo-axis in transverse to strain 
contact configuration in stress-free (blue line) and stressed 
(partially detwinned, yellow line) states. Inset shows zoom of 
the transition area in normalized resistance plot in compari- 
son with resistance of sample D in the detwinned state. 



anisotropy drops to 1 right above the transition, reflect- 
ing its 1st order nature, in Bal22 an anisotropy is de- 
tectable at least up to 30 K above the transition, a man- 
ifestation of strong nematic fluctuations above Tto • 
This is consistent with the phase transition in Bal22 be- 
ing second order or very close to such, which, in term, 
was argued to be related to a harder lattice and weaker 
magnetoelastic coupling in Bal22. 



X-ray characterization of detwinned state of 
BaFe 2 As2 



FIG. 5: (Color online) Top panel. Image of the sample for 
resistivity measurements in transverse to the strain (horizon- 
tal) direction. The current contacts are made along the whole 
length of the top and bottom edges of the sample, right and 
left sample edges are covered with with epoxy glue to provide 
homogeneous deformation of the sample and wires. Straining 
force is applied on horseshoe, similar to Fig. lc. Potential 
contacts are made small to minimize the effect on the strain 
distribution. Middle and bottom panels show area close to 
potential contacts in stress-free and stressed states. 



that this similarity holds despite the fact that the resid- 
ual resistivity in Bal22 is substantially higher — the 
anisotropy of po is also higher. As explained below, this 
is consistent with the idea that the transport anisotropy 
reflects the anisotropy of the carrier velocity. 

Another important observation is that in Cal22 the 



The upper panels of Fig. [7J show the splitting of 
the tetragonal (220)t reflection into the orthorhombic 
(400)oi and (040)o2 reflections produced by domains Ol 
and 02, with different orientations, below the structural 
transition in a sample area that was nearly strain free. 
Illuminating the sample with 200/im diameter x-ray spot 
at positions between strain-applying contacts yields the 
pattern shown in the lower panels of Fig. [7] The intensity 
in the (400)oi reflection is strongly enhanced indicating 
that the domains with their orthorhombic a axis (a > b) 
along the direction of the applied strain are more popu- 
lated. The full penetration of the high-energy x-ray beam 
through the sample allows a quantitative analysis of the 
bulk domain population with an error of approximately 
1 % by using the relative intensities in all orthorhombic 
(400) and (040) reflections. 

Scanning the x-ray beam across the sample allows a 
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FIG. 8: Spatial distribution of the population of Ol domains 
in strained BaFe2As2. The population was determined from 
the relative intensities of the (400)oi and (040)o2 peaks, as 
shown in Fig. [7] The area between the contacts, through 
which the strain was applied, was probed in 0.2 mm vertical 
and horizontal steps, determined by the :r-ray beam size. The 
number and color in each pixel correspond to the percent- 
value for the domain Ol population. The pixel in the upper 
right corner corresponds to an area under the contact. 



FIG. 7: Sections of reciprocal (HK0) planes around the 
position of the tetragonal (220) t reflection recorded from 
the BaFe2As2 sample by high-energy x-ray diffraction as de- 
scribed in the text. In the orthorhombic state at T = 25 K, 
the (400)oi and (040)o2 reflections are related to different do- 
mains Ol and 02, respectively. Analysis of the intensities in 
pattern recorded at different positions on the sample, below a 
contact in the upper panels and between the contacts in the 
lower panels, yield the indicated %- values for the domain Ol 
population bolstered by the applied strain. 



spatially resolved characterization of the domain pop- 
ulation as demonstrated in Fig|8j Here we show that 
the sample is nearly detwinned between the contacts, 
whereas the domain population comes closer to 50% (the 
value expected for a random distribution of twins) at 
positions below the contacts, as shown in the right-top 
pixel in Fig [8] This confirms our conclusion, Fig. [2l that 
contacts act as anchors for the deformation. 



Theoretical analysis 

To analyze these results let us write the resistivity in 
Boltzmann (relaxation time) approximation, keeping in 
mind that the relaxation times for holes and for elec- 
trons may be different. Despite the fact that the mag- 
netic structure of A122 compounds is highly anisotropic, 
we would assume that the relaxation times is isotropic, 



as we will discuss below (in particular, it is important 
to keep in mind that the Fermi surfaces, both experi- 
mental and calculated, are much smaller in size that the 
antiferromagnctic vector). Thus, we write 



VP = VPh + VPe 



(1) 



Pi = (47r/wp j)(l/To,i + l/r p h,i + l/r ee ,i + l/r mofl) i), 

(2) 

where i stands for hole or electrons. The prefactor, wj, 
is the only quantity directly affected by twinning (except 
for a possible scattering by domain walls in tq that in 
principle would go away with detwinning; however, the 
twin domain walls are usually weak scatterers). In the 
classical Boltzmann theory it is presumed that lo^ and To 
are temperature- independent, which allows one to sepa- 
rate p additively into two terms, a temperature indepen- 
dent residual resistivity, and the temperature dependent 
scattering by thermal excitations (phonons, magnons, 
electrons). Unfortunately, this is not the case here. The 
Hall effect studies found that upon cooling in the AFM 
phase, the carrier concentration drops drastically, as re- 
flected by the nearly an order of magnitude reduction 
in the plasma frequency, tic, and the density of states 
also decreases substantially (the latter effect makes To 
temperature dependent, but this is not important for us 
now). 

The standard approximation implies an isotropic re- 
laxation time, even when the bosons providing inelas- 
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tic scattering arc anisotropic. Our first task is to ver- 
ify whether this remains a good approximation in our 
case. Given the large magnetic anisotropy and large dif- 
ference in the ferro- and antiferromagnetic correlation 
length above Tto [34 - 36 1 one may think that scattering 
by the spin fluctuations may be anisotropic as well. To 
analyze that, let us recall the classical theory of the in- 
elastic transport scattering in solids. The scattering rate 
for the current flowing in a given direction i is defined by 
the Fermi golden rule (see, e.g., Ref. |37J, Eq. 37]): 



k 



kk' 



where A is an isotropic factor, the denominator is pro- 
portional to the plasma frequency squared, is the 
matrix element of the electron-magnon scattering by a 
magnon with the wave vector q = k — k' and energy 
hi>j(c\). and the cross-terms involving VkVk> were ne- 
glected. The standard argument goes like this: even 
though the bosons (magnons, in our case), may be much 
softer, and more easily excited by temperature, in one 
g-dircction than in the other, the fact that the scattering 
processes are integrated twice over the Fermi surface ef- 
fectively averages out any spectral anisotropy. Moreover, 



in fact direct neutron measurements [35l . l36j ] show that 
the magnon dispersions along x and along y differ by less 
than 15%. Furthermore, small size of the Fermi surfaces 
in the magnetic state (note that the Hall concentration 
drops quite rapidly below Tto , indicating that the Fermi 
surfaces start to shrink immediately below the transition 
[38j |) only lets the magnons with small q participate in 
scattering, and at small qs the anisotropy of the magnon 
spectra is particularly low. Thus a strong anisotropy in 
r is unlikely, although cannot be excluded. 

Let us now analyze the most anisotropic term in Eq. 
[H the plasma frequencies. These are intimately related 
with the electronic structure near the Fermi level, so we 
show in Fig. [5] the calculated Fermi surfaces for both 
compounds in their calculated ground states. Note that 
while the calculated ground state moment is larger than 
the experimental, the calculated number of carriers (the 
Fermi surface volume) is also larger that the experimental 
one H, 39|. Therefore it is hard to say what would be 



a better calculation to use for the transport properties: 
GGA gives more correct (but still too large) volumes of 
the Fermi surfaces, while LDA gives more correct (but 
still too large) magnetic moments. Besides, as we will 
show next, the calculated anisotropy is very sensitive to 
these details, primarily because of the small size and the 
complicated shape of the Fermi surfaces. That is to say, 
the computational results should be taken with a large 
grain of salt, and even if they provide in some aspects 
reasonable quantitative agreement with the experiment, 
this is likely fortuitous and one should not be tempted to 




FIG. 9: (color online) Fermi surfaces for Cal22 (left panels) 
and Bal22 (right panels), obtained in the LDA band struc- 
ture calculation (low magnetic moment). One reciprocal lat- 
tice unit cell is shown, with the T points located in the eight 
corners. The top panels represent a general view (directions 
of the axes as shown), the next row represents the view from 
the top, along the k z direction (k x -k y plane), the third row 
represents the front view, along the k y direction (k x -k z plane), 
the last row presents the side view, along the x direction (k y - 
k z plane). The light shaded (cyan) surfaces are hole pockets, 
the dark (blue) one are the electron pockets. 



use the calculated numbers any more than a qualitative 
guide. 

Nevertheless, one can readily make some qualitative 
observations: first, both electron and hole Fermi sur- 
faces strongly violate the tetragonal symmetry; second, 
the two com pou nds differ substantially (this has been no- 
ticed before [40[ and ascribed to the direct As- As overlap 
across the Ca plane). One can naively think that the 
calculated in-plane transport anisotropy should be very 
large, yet it appears not to be the case (Table 1). Partic- 
ularly, in the antiferromagnetic LDA ground state Cal22 
has an anisotropy a xx /a yy of 1.10 (+10%), and in the 
GGA ground state 0.97 (-3%). In Bal22 the correspond- 
ing numbers are +14% and +31%. 

Interestingly, the reason for such low anisotropy is dif- 
ferent for the two compounds. As one can see from the 
Table 1, in Cal22 the holes and the electrons separately 
have noticeable anisotropics of the opposite sign: —42% 
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TABLE I: Calculated total and partial (hole, electrons) 
squared plasma frequencies in eV 2 . The x axis is parallel 
to the crystallographic a axis and perpendicular to the fer- 
romagnetic chain direction. The second column indicates the 
calculated magnetic moment in [lb /Fe. 

M u? vx uly iol x {h) toly(h) Upjfi) ul v (e) 



Cal22 (GGA) 
Cal22 (LDA) 
Bal22 (GGA) 
Bal22 (LDA) 



1.85 0.36 0.37 0.14 0.24 0.22 0.14 

1.61 0.46 0.42 0.11 0.24 0.35 0.18 

1.96 0.42 0.32 0.13 0.13 0.29 0.19 

1.64 0.41 0.36 0.20 0.17 0.20 0.20 



and +57%, respectively, in GGA, and -46% and +94% 
in LDA. The net anisotropy is thus small owing to the 
large cancellation of these two terms. This cancellation 
holds only if the relaxation times for holes and electrons 
are similar. The validity of this approximation in Cal22 
is unknown, but in Bal22 it is known to be highly ques- 
tionable. Based on the Hall effect data, it was argued 



39J that the electrical transport in Bal22 at low tem- 
peratures is dominated by electrons. 

On the other hand, in Bal22 both the holes and elec- 
trons separately show rather weak anisotropy of the same 
sign, especially in LDA: 18% for the holes and 0% for 
the electrons. In GGA, these numbers arc and +52%. 
The reason is that while the Fermi surface pockets in 
Bal22 are aligned along x, violating the tetragonal sym- 
metry, each of these pockets separately (and as opposed 
to Cal22) shows relatively little anisotropy. 

While these numbers are not quantitatively reliable, 
they show two qualitative trends observed in the ex- 
periment: firstly, despite the fact that the Fermi sur- 
faces themselves are strongly anisotropic, the transport 
anisotropy is weak in Cal22 and moderate in Bal22. Sec- 
ondly, despite the fact that the magnetic structure sug- 
gests much higher conductivity along y (i.e, along the fer- 
romagnetic chains) than along x, the opposite holds both 
in the calculations and in the experiment. 

The temperature dependence is more intriguing. Ob- 
viously, if one could plot the anisotropy as a function of 
the long-range staggered magnetization, in has to be fi- 
nite when the magnetization is maximal, at T = 0, and 
zero when the magnetization is vanishingly small. Yet 
in the experiment the anisotropy is growing with tem- 
perature, while the moment decreases. This means that 
right below Tto the magnetic moments are already large, 
even though incompletely ordered (as oppose to what the 
fully spin-Peierls model 41, 42| would imply) and so is 
the transport anisotropy. An example of such model is 
the dynamic domain scenario [23} that suggests that the 
long-range order is destroyed primarily by fluctuating do- 
main walls that affect the average long-range magnetiza- 
tion, but do not affect electronic anisotropy (see Fig. la 
in Ref.[23t). 

Within this scenario we can understand why the 



anisotropy experiences a nearly finite jump below TtOj 
however, one needs to explain why the anisotropy ac- 
tually decreases with cooling. One possible explanation 
is that the ratio of the plasma frequencies nonmonotoni- 
cally depends on the magnitude of the magnetic moment. 
In this case anisotropy will be correlated with the mo- 
ment for very small (experimentally inaccessible) mag- 
netization, but anticorrelated with the moment at larger 
moments. 

Alternatively, one may think that if the hole and 
the electron anisotropics have opposite signs, the net 
anisotropy is defined by their incomplete cancelation. In 
this case the observed temperature dependence is de- 
fined by the temperature dependence of the Th/r e ratio, 
which, of course, can have virtually any temperature de- 
pendence. This latter scenario is corroborated by the 
calculations in Cal22, but not in Bal22. One has to re- 
member, as we mentioned before, that these calculations 
are not to be trusted on a quantitative level. 



CONCLUSIONS 

In conclusion, we performed reversible mechanical de- 
twinning of single crystal CaFe2As2 and BaFe2As2. The 
single domain state was confirmed by direct polarized 
light imaging in both compounds and with high-energy 
x-rays in BaFe2As2. In both materials we find that the 
resistive anisotropy is largest at Tto with p ao /pbo ~ 1-2 
in Cal22 and wl.5 in Bal22. For Cal22 this anisotropy 
only exists below Tto and diminishes upon further cool- 
ing, reaching about 1.05 at T ~ 50 K and remaining 
around that value at lower temperatures. For Bal22 
the anisotropy exists both below and above Tto , reflect- 
ing nematic fluctuations above Tto an d the second-order 
character of the phase transition there. The temperature 
dependence of the anisotropy is weaker in Bal22, the 
low-temperature value being comparable with the high- 
tcmpcrature one, 1.4 vs. 1.5. 

These results are rather counterintuitive in several as- 
pects. First, if one adapts a simple Jahn- Teller-like 
orbital ordering picture [4314461 ] . one expects that con- 
ductivity along the ferromagnetic "metallic" chains (y) 
should be higher than along the antiferromagnetic "in- 
sulating" direction (x). The opposite is true. Second, 
given the drastic anisotropy of the band structure and 
the Fermi surfaces (including the fully broken symmetry 
between the xz and yz orbitals), as confirmed by both 
first principles calculations and the Hall experiments, one 
expects a dramatic anisotropy, while in fact we observe 
weak to moderate anisotropy, not atypical for many met- 
als. Third, within the experimental accuracy, the max- 
imal anisotropy is right below the transition, and not 
when the order parameter (magnetic moment or the or- 
thorhombicity) is maximal. Finally, despite the fact that 
the degree of orthorhombicity and the magnetic order pa- 



10 



ramctcr monotonically decrease with temperature, the 
transport anisotropy actually increases. 

The first two findings are fully corroborated by the first 
principle calculations. These predict the correct (coun- 
terintuitive) sign of the anisotropy, and a relatively mod- 
erate anisotropy magnitude. The third might reflect the 
existence of dynamic domains (similar to the nematic 
fluctuations) right below T TO in Cal22, and both below 
and above in Bal22. In that case, the antiphase domain 
walls, as described in (23|, break the long range order 
without destroying the transport anisotropy. 

As to the temperature dependence, we have no prefer- 
able interpretation of this effect; the most plausible cause 
seems to be the temperature dependence of the relaxation 
rates, which are known to change by more than an order 
of magnitude between T = and T = Xto- 

While this manuscript was prepared for submission, 
a preprint appeared reporting a partial mechanical de- 
twinning using uniaxial stress in Ba(Fei_ :r Co^ As 2 [47 1. 
The reported resistivity data for the parent BaFe2As2 are 
similar to those found in our study. 
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